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ABSTRACT 

The domain wall between coexisting chirally symmetric and broken— symmetry 
regions is studied in a saddle— point approximation to the effective three— flavour c— 

model. In the chiral limit the surface tension varies in the range [(40 50)MeV} 3 . 

The width of the domain wall is estimated to be ~ 4.5 fin. 
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1. Recently, considerable effort has been invested into the determination of the 
surface tension between the low- and the high-temperature phases of quantum 
chromodynamics at the coexistence point [1,2]. These investigations have made 
use of the quenched lattice version of QCD. Already at modest temperature- 
extension of the lattice (iVj = 4) considerable difficulties were met in extracting a 
reliable signal for the presence of the interface. In the full theory, the simulation 
of the phase transition with 4 light flavours is well under control [3]. In the 
phenomenologically more interesting Nf =2,3 cases the order of the transition is 
not settled yet [4,5]. 

On the basis of these remarks the renewal of interest in the effective low 
energy theories of QCD seems to be justified in the present context too. Since 
some exploratory investigations indicate that near the phase transition the high 
temperature phase of QCD might also be described in terms of colourless, massive 
excitations [6], there is a chance that the phase transition itself can be understood 
fully in this framework. 

In this Letter we present a semiquantitative analysis of the interface between 
coexisting chirally symmetric and broken symmetry regions in a three-flavour 
linear <r-model. 

The most successful effective meson model is based on the non-linear realisa- 
tion of the broken-symmetry vacuum [7,8]. In particular, based on the observation 
that the finite-temperature partition function of QCD is determined at low T by 
the Goldstone sector, an expansion has been derived for the quark condensate, 
to 0(T 6 ) [9]. A smooth extrapolation of its decreasing tendency towards 
higher T indicates that the restoration of chiral symmetry would take place at 
« 170 MeV for Nf = 2. In view of the possible first-order nature of the transition 
this is rather an upper bound on T c . As such, it is in agreement with the results of 
MC— simulations of the full QCD with two [10] and four [3] light flavours (in the lat- 
ter case the proportionality of T c to N j is to be taken into account). However, 

a fully conventional approach to the chiral phase transition within the non-linear 
model is prevented by the absence of an explicit order parameter field in this the- 
ory. Actually, this field is traded for the coupling F wt whose temperature- variation 
has been also studied in the literature [9,11]. 

An interesting variant of this low-energy description incorporates also a glue- 
ball field, fulfilling the trace anomaly already at the tree level [12,13,14]. The 
resulting idea of a unique mechanism for the finite temperature variation of the 
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chiral and gluon condensates has been investigated recently in Ref. [15], 

It is in the linear <r-model where the scalar singlet order parameter field 
appears explicitly. Although the existence of the scalar nonet is a notoriously 
disputed subject, recent analyses claim, that by taking into account the existence 
of K-K molecules, a nonet with the appropriate features can be proposed uniquely 
[16]. The linear model at the tree level represents equally well the low-energy 
meson sector as the non-linear version [17,18,19]. However, it has been explicitly 
shown that renormalisation is not the correct way to improve its predictions [20], 
Therefore, it is a reasonable strategy [21] to estimate the modification of the tree- 
level potential at finite temperature by taking into account exclusively the thermal 
fluctuations. 

This “classical” interpretation of the effective model has been implemented in 
Ref. [21] by simply omitting the T-independent (divergent) part of the one-loop 
finite temperature quantum correction of the effective potential. This certainly 
cannot be considered a consistent approach. Instead, we emphasize that the cr- 
model provides a classical field-theoretical model for the ground state of QCD, 
which can be investigated at finite temperature in the framework of classical sta- 
tistical physics. We remark that recently the role played in the deconfinement 
by skyrmions has been investigated in the SU( 2) x 5C/(2) non-linear model with 
the same philosophy [22]. Before proceeding to the construction of the chiral 
order parameter profile, we shall calculate the constrained free energy [23,24] of 
the tree-level parametrized linear <r-model as for a classical field theory at finite 
temperature. A first-order transition will be detected for Nf = 3 in agreement 
with the results of the general renormalisation group analyses [25,26]. Once the 
functional form of the free energy constrained by the constant non-zero value 
of the <r-field is calculated at T c numerically, the construction and the analysis 
of the chiral order parameter profile goes along the same lines as earlier for the 
Polyakov-loop profile [27]. 

The present approach is in full analogy with the treatment of magnetism with 
the help of classical spin-models (Ising, Heisenberg etc.), where the quantal nature 
of the magnetic phenomena is reflected only in the actual value of the nearest 
neighbour, etc. couplings calculated at zero temperature. The applicability of 
such models to real magnetic crystals depends on the success or failure of their 
predictions. The same statement applies to the content of the present note. 
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2. The Hamiltonian of the linear cr-model is given by 

H[M] ^ x tr ■*" 

H P [M ] = H£ tr (VM)(VM^ ) - y tr MAft + g{ detAf + det Aft) 

+/j( tr AfAft) 2 + f 2 tr (AfAft) 2 , (1) 

where the complex matrix M represents the scalar (<x) and the pseudoscalar (</>) 
nonets: 

-i 8 

M(x) = —p + i<f>i{ x )) A/, tr A|Ajt = 26 ifc . (2) 

Isb0 

The matrix II is formed from the conjugate momenta of M. In our concrete 
calculations we shall make use of the parametrisation proposed by Goldberg [21], 
which leads to an ££/(3)— symmetric spontanous breaking of the £1/(3) X £1/(3) 
chiral symmetry: 

fi 0 = 244MeV, g = -1811MeV, fi = 1.70, f 2 = 11.91. (3) 


For the construction of the chiral profile the constrained free energy, corre- 
sponding to a non-zero fixed value of <To(k = 0), the k = 0 Fourier component of 
<ro(z), is the appropriate quantity to evaluate. After integrating over the conjugate 
momenta the partition function of the classical medium defined by Eq. (1) and 
Eq. (3) is given by 

Z[cr 0 ] = exp{-(3VU((T 0 )} 


J VM6{a 0 { k = 0 ) 


8 8 

- o-o ) n = °)) n^W k = °)) ex P {-PUp}, 

1=1 1=0 


( 4 ) 


where (3 = T 1 . 

The integration over M in a finite volume V is done by separating the k = 0 
mode explicitly: 


M(x) = —y=(T 0 I + m(x), m(k = 0) = 0. (5) 

\/3 

Reshuffling Eq. (1) accordingly, one proceeds to an approximate evaluation of the 
m— integral, which would be exact in an infinite component model (actually, here 
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one has 18 components). First, an auxiliary 3x3 matrix field E(x) is introduced, 
which formally transforms the quartic piece in the exponent of Eq. (4) into one 
which is quadratic in rn . Next, one assumes that the E— integration is dominated 
by an 5I7(3)-symmetric saddle-point: 

E(x) = i I. (6) 

Third, one performs the m-integration in the Gaussian approximation, which 
yields the leading low-temperature correction to the zero-temperature free en- 
ergy. (The perturbation theory with respect to the cubic piece of Tip corresponds 
to a loop expansion, where the role of h is played by /? _1 = T.) 

The detailed presentation of the above procedure can now be found in text- 
books [28], therefore we directly give the infinite-volume result for the one-loop 
correction of the constrained free energy, evaluated at the saddle point: 

yAU.' idl . =Tj J2 »(«) + *«)) ' ( 7 ) 

In this equation the index Q runs over the scalar singlet (SO), the scalar octet (S8), 
the pseudoscalar singlet (PO), and the pseudoscalar octet (P8) meson-multiplets, 
g(Q) gives the corresponding multiplicities, X{Q) = s + fi 2 + fx 2 (Q) acts as a kind 
of effective squared mass, while /z 2 (Q) at <7 0 = <r m j n (T = 0) yields the meson 
masses in the chiral limit: 

/i 2 (S0) = — /z 2 4- — + 4 u 2 (3 fi + fo), 

^ 2 (58) = — /z 2 — — ^=«7o + (/i + fi ), 

/i 2 (P0) = -n\ - ^j=<T 0 + ^(3/j + /2), 

/x 2 (P8) = ~nl + + /2), 

One has to emphasize that the original integrand of Eq. (4) is exponentially small 
for large field values, therefore the saddle-point approximation loses its sense when 
any single X(Q) becomes negative. In the symmetric phase the particle masses 
are defined at <To = 0, where i has the meaning of the common mass squaxe of 
the parity partners. At zero temperature i = — /z 2 (see below), which implies by 


5 (50) = 1, 
5(58) = 8, 
S(P0) = 1, 
U(P8) = 8. 


( 8 ) 
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Eq. (3) & min = 115 MeV leading to a Sl7(3)-symmetric spectrum (P8: 0, P0: 850 
MeV, S8: 950 MeV, SO: 600 MeV). Thus, the saddle-point approach holds the 
promise for a unified description of both phases. 

The saddle-point value of the free energy is found by minimising its expression 
as a function of J for fixed ctq : 


9t^taddle{^ ) ^o] q 

ds 


( 9 ) 


The formal procedure outlined above is given sense by applying a cut-off A 
on the k-integral in Eq. (7). This is physically natural, since the effective theory 
ceases to represent QCD, even approximately, when high-energy fluctuations start 
to play an important role. For the numerical calculations it is convenient to use 
scaled (dimensionless) quantities: 


T = 


T 

A’ 





( 10 ) 


The cut-off integral of Eq. (7) was calculated analytically: 


[-,«r, A ] = A4 {~ 8 (3/ZV/,) ^ 2 + ~ I* 1 *' + + (/l + ^ 

+ + 2 X(Q) - 2X(Q) 3 / 3 arctan X~ 1/2 (Q)] J. 

Q (11) 
The explicit expression of Eq. (9) can be obtained by taking the partial derivative 
of Eq. (11) with respect to s . 

Eq. (9) has been solved for several values of A in the range A G (250, 2000)MeV . 
For each A, r was varied and s(A, r, <r) found. In the range A G (500, « 1000) MeV 
for some <r-values no saddle-point fulfilling the X > 0 condition could be found. 
In these points the proposed approximation scheme fails. Based on the T c -range 
of lattice-calculations discussed in the introduction, our strategy was to choose 
values of A tuning the critical temperature over the interval T c £ (100,200) MeV. 

It turns out that the interesting interval in this sense is A < 400 MeV, where the 
saddle-point approximation works well. _ 

Evaluating Eq. (11) at the saddle point, at very low t its double— well nature 
inherited from the classical potential is preserved. When increasing t, the position 
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of the non-trivial minimum shifts considerably towards the symmetric minimum 
at the origin. The degeneracy of the corresponding potential values occurs before 
their locations coincide. In the whole parameter region, where the saddle— point 
approximation makes sense, a clean first-order transition was observed. It is worth 
mentioning that on both sides of the transition temperature we found physical 
values for the saddle point ( X > 0), which can be contrasted with the case of 
continous phase transitions, discussed in [28], where the critical point is identified 
from the symmetric phase as the end of the validity of the saddle-point dominance. 

In Fig. 1 we illustrate the typical evolution of the potential with the tem- 
perature for A = 300MeV. Other values of the cut-off lead only to quantitative 
changes. In Fig. 2 the smooth variation of T c with A is displayed. Although 
T c < A in the whole range, the relevant A values are lower than what is intuitively 
expected. We remark that T c < lOMeV for A > lGeV. 

3. When the critical temperature is fixed by choosing A, no free par am eter remains 
in the theory. The complete characterisation of the chiral interface will be given 
in physical units. The kink equation 

<fi ! da,, { l 

is discretised using a lattice constant (a) proportional to A -1 


2J iQf , CL — j 71 — 1, 2 j • •• 

71 A 


(13) 


The discretised equation for the dimensionless quantities was solved in the form 


.(( + !) = 2.(1) -,(,-1)+-^!^. (14) 

The kink solution of the A = 300 MeV case suggests very slight n— dependence, 
even for n = 1,2 (Fig. 3). In addition, we have calculated the surface tension, 


= E [\ n M l ) - - !)) 2 + Aj W..«ie(<r)] (15) 

which also shows very slight n-dependence. In Table 1 we give the surface tension 
(a), the width of the interface (<f) and the value of the non-trivial <r m j n at T c for 
a few T c -values from the relevant range. The most important observation is that 
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none of them seem to be particularly sensitive to T c . The results summarized in 
the Abstract are rather robust. 

The range of a found by us is almost one order of magnitude smaller than the 
bag-model estimate of Ref.[29]. Although Campbell et al have given some support 
to that result in their effective model calculation, our thickness estimate sheds 
some doubt on the validity of the thin— wall approximation used in [15]. This wall 
thickness of 4—5 fm is a somewhat intriguing result of the present study. It would 
mean that in heavy ion collisions, in a large part of the fireball only a decrease 
of the quark— condensate could be achieved, as the tail of the surrounding low- 
temperature vacuum enters fairly deep into the high energy— density region. Such 
a situation might motivate the study of particle emission from states intermediate 
between the fully broken and the symmetric phases at T c . 

Our additional numerical investigation showed that if the mass of the scalar 
singlet is increased by up to 950 MeV (the other parameters are kept fixed) no 
quantitative change occurs in the characteristics of the interface. However, if the 
mass of the whole scalar sector is larger (of the order of 1.5 GeV) the surface 
tension would grow to the value conjectured in Refs. [32], and its width would 
decrease to 2 frn. 

4. In this investigation we have attempted a semi-quantitative characterisation 
of the chiral interface at the QCD phase transition point. The linear cr-model 
was just a choice, where the temperature dependent corrections to the free energy 
density of the er— field could be estimated. The model certainly has to be extended 
to take into account further light degrees of freedom with significant coupling to a . 
A first step in this direction could be a model incorporating an additional glueball- 
field [14]. We note, that with careful discretisation, which takes into account the 
essential role played by the cut-off, one can use non— perturbative lattice methods 
for the study of the thermodynamics of low-energy effective theories as well. 
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Table 1 


The transition temperature and the characteristic data of the chiral order 
parameter profile in the chiral limit of the linear <r-model as a function of the 
cut-off. 


A [MeV] 

T c [MeV] 

a 1 / 3 [MeV] 

d [fm] 

Cmin [MeV] 

230 

225.77 

38.9 

4.5 

48.5 

240 

199.14 

39.7 

4.5 

50.2 

250 

176.66 

40.6 

4.5 

51.8 

260 

157.55 

41.3 

4.5 

53.0 

270 

141.20 

42.2 

4.5 

54.3 

280 

127.14 

43.1 

4.5 

55.7 

290 

114.97 

43.9 

4.5 

56.8 

300 

104.38 

44.6 

4.4 

57.9 

310 

95.13 

45.4 

4.4 

58.9 

320 

87.01 

46.1 

4.4 

59.6 


Figure captions 


Figure 1: The variation of 10 3 x the scaled constrained free energy density vs. the 
scaled c-field with the temperature. The curves, starting from below, 
correspond to T = 100 MeV, 102 MeV, T c , 106 MeV and 108 MeV, 
respectively. 

Figure 2: The dependence of the critical temperature on the cut-off of the thermal 
fluctuations. Both T c and A are measured in MeV. 

Figure 3: The chiral order parameter profile displayed on a lattice with lattice con- 
stant a=l/3A. The different curves correspond to appropriately scaled- 
up kinks evaluated originally on lattices with a=l/A (solid), a=l/2A 
(dashed) and a=l/3A (dotted), respectively. The value of the <r-field is 
measured in MeV. 
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